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Abstract 



We present a measurement of the transverse momentum distribution of Z 
bosons produced in pp collisions at -y/s = 1.8 TeV using data collected by the 
D0 experiment at the Fermilab Tevatron Collider during 1994-1996. We find 
good agreement between our data and a current resummation calculation. We 
also use our data to extract values of the non-perturbative parameters for a 
particular version of the resummation formalism, obtaining significantly more 
precise values than previous determinations. 
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We report a new measurement of the differential cross section with respect to trans- 
verse momentum (da/dpT) of the Z boson in the dielectron channel with statistics and pre- 
cision greatly improved beyond previous measurements The measurement of da/dp-r 
of the Z boson provides a sensitive test of QCD at high-Q^. At small transverse momentum 
(Pt), where the cross section is highest, uncertainties in the phenomenology of vector boson 
production have contributed significantly to the uncertainty in the mass of the W boson. 
Due to its similar production characteristics and the fact that the decay electrons can be 
very well-measured, the Z provides a good laboratory for evaluating the phenomenology of 
vector boson production. 

In the parton model, Z bosons are produced in collisions of qq constituents of the proton 
and antiproton. The fact that observed Z bosons have finite px can be attributed to gluon 
radiation from the colliding partons prior to their annihilation. In standard perturbative 
QCD (pQCD), the cross section for Z boson production is calculated by expanding in pow- 
ers of the strong coupling constant, as- This procedure works well when p'^ ~ with 
Q = Mz- However, when pr ^ Q, correction terms that are proportional to a<( ln((5^/p^) 
become significant, and the cross section diverges at small px- This difficulty is surmounted 
by reordering the perturbative series through a technique called resummation [p|-[T3[|. Al- 



though this technique extends the applicability of pQCD to lower values of px, a more 
fundamental barrier is encountered when px approaches Aqcd- In this region, as becomes 
large and the perturbative calculation is no longer valid. In order to account for the non- 
perturbative contribution, a phenomenological form factor must be invoked, which contains 
several parameters that must be tuned to data [3,P^,P^ ■ 



The resummation may be carried out in impact-parameter (6) space via a Fourier trans- 
form, or in transverse momentum space. Both formalisms require a non-perturbative func- 
tion to describe the \ow-px region beyond some cut-off value bmax or pxiim and they merge 
to the fixed-order perturbation theory at px ~ Q- The current state-of-the-art for the 
6-space formalism resums terms to next-to-next-to-next-to-leading-log and includes fixed- 
order terms to 0{a'^) [[l^]. Similarly, the p^-space formalism resums terms to next-to-next- 
to-leading-log and includes fixed-order terms to 0{as) [jlSj . 

In the 6-space formalism, the resummed cross section is modified at large b (above bmax) 
by exp(— S'np(&, Q^)). The form factor S'np(&, Q^) has a general renormalization group in- 
variant form, but requires a specific choice of parameterization when making predictions. A 
possible choice, suggested by Ladinsky and Yuan [|ll|, is 



gib"" + g2bHn{^J + g,g3b\n{100xiXj), (1) 

where Xi and Xj are the fractions of incident hadron momenta carried by the colliding 
partons and gi are the non-perturbative parameters. An earlier parameterization by Davies, 
Webber, and Stirling |g| corresponds to the above with (73 = 0. For measurements at the 
Fermilab Tevatron at = M|, the calculation is most sensitive to the value of g2 and quite 
insensitive to the value of gs. 

In the px-space formalism, the resummed cross section is modified at low-px (below pxnm) 
by multiplying the cross section by -Fnp(pt)- In this case, the form of the non-perturbative 
function is not constrained by renormalization group invariance. The choice suggested by 
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Fnp(pt) = 1 - e-'^'PT 



(2) 



where a is a non-perturbative parameter. 

Previously published measurements of the differential cross section for Z boson produc- 
tion have been limited primarily by statistics (candidate samples of a few hundred events). 
This measurement is based on a sample of 6407 Z — » e^e~ events, corresponding to an 
integrated luminosity of ~ 111 pb~^, collected with the D0 detector in 1994-1996. A 



recent measurement by the CDF Collaboration has a similar number of events |]15 



Electrons are detected in the uranium/liquid- argon calorimeter with a fractional energy 



resolution of ~ 15%/ y -E'(GeV). The calorimeter has a transverse granularity at the electron 
shower maximum of Ar] x A0 = 0.05 x 0.05, where rj is the pseudorapidity and is the 
azimuthal angle. The two electron candidates in the event with the highest transverse 
energy (-Er), both having Et> 25 GeV, are used to reconstruct the Z boson candidate. One 
electron is required to be in the central region, |?7det| < 1-1, and the second electron may 
be either in the central or in the forward region, 1.5 < |?7det| < 2.5, where rj^et refers to the 
value of 7] obtained by assuming that the shower originates from the center of the detector. 
Offline, both electrons are required to be isolated and to satisfy cluster-shape requirements. 
Additionally, at least one of the electrons is required to have a matching track in the drift 
chamber system that points to the reconstructed calorimeter cluster. 

Both the acceptance and the theory predictions modified by the D0 detector resolution 
are calculated using a simulation technique originally developed for measuring the mass of 



the W boson |]T6[, with minor modifications required by changes in selection criteria. The 
four-momentum of the Z boson is obtained by generating the mass of the Z according to an 
energy-dependent Breit-Wigner lineshape. The pt and rapidity of the Z boson are chosen 



randomly from two-dimensional grids created using the computer program LEGACY [12 



which calculates the Z boson cross section for a given p-r, rapidity, and mass of the Z 
boson. The positions and energies of the electrons are smeared according to the measured 
resolutions, and corrected for offsets in energy scale caused by the underlying event and 
recoil particles that overlap the calorimeter towers. Underlying events are modeled using 
data from random inelastic pp collisions of the same luminosity profile as the Z boson sample. 
The electron energy and angular resolutions are tuned to reproduce the observed width of 
the mass distribution at the Z-boson resonance and the difference between the reconstructed 
vertex positions of the electrons. 

We determine the shape of the efficiency of the event selection criteria as a function 
of Pt using Z e^e~ events generated with herwig [|I7|, smeared with the D0 detector 
resolutions, and overlaid on randomly selected zero bias pp collisions. This simulation models 
the effects of the underlying event and jet activity on the selection of the electrons. The 



absolute efficiency is obtained from Z — > e~^e~ data [|T^]. The values of the efficiency times 
acceptance range from 26-37% for px below 200 GeV and is 53% for px above 200 GeV. 

The primary background arises from multiple-jet production from QCD processes in 
which two jets pass the electron selection criteria. We use several D0 data sets for esti- 
mating this background — direct-7 events, dijet events, and dielectron events in which both 
electrons fail quality criteria — all of which have very similar kinematic characteristics [|l|. 
The level of the multijet background is determined by fitting the ee invariant mass in the 
range 60 < Mgg < 120 GeV to a linear combination of Monte Carlo Z —* e^e~ signal 
events (using pythia |T^) and background (from direct-7 events). We assign a systematic 
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FIG. 1. The differential cross section as a function of px compared to the resummation calcu- 
lation with three different published parameterizations of the non-perturbative region and to the 
fixed-order calculation. Also shown are the fractional differences between the data and each of the 
resummation predictions. 



uncertainty to this measurement by varying the choice of mass window used in the fit, and 
by changing the background sample among those mentioned above. We estimate the total 
multijet background level to be (4.4±0.9)%. The direct-7 sample is used to parameterize the 
shape of the background distribution as a function of pt ■ Backgrounds from other sources, 
such as Z — >• r^r~, tt, and diboson production, are negligible. 

We use the data corrected for background, acceptance, and efficiency, to determine the 
best value of the non-perturbative parameter, g2, given our data. In the fit, we fix gi and 
(?3 to the values obtained in |ll[] and vary the value of g2- We use the CTEQ4M pdf. The 
prediction is smeared with the known detector resolutions, and the result fitted to our data. 
The resulting distribution as a function of g2 is well-behaved and parabolic, yielding a 
value of g2 = 0.59±0.06 GeV^, considerably more precise than previous determinations. For 
completeness, we also fit the individual values of gi and gs, with the other two parameters 



fixed to their published values We obtain gi = 0.09 ± 0.03 GeV and g^ = —1.1 ± 0.6 



GeV~^. Both results are consistent with the values of Ref. [IT 

To determine the true da/dpT, we correct the measured cross section for effects of detector 
smearing, using the ratio of generated to resolution-smeared ansatz pj- distributions. We 
use the calculation from legacy as our ansatz function, with the g2 determined from our 
fit. The largest smearing correction occurs at low-pr, where smearing causes the largest 
fractional change in pt and where the kinematic boundary at pt= produces non-Gaussian 
smearing. The correction is 18.5% in the first bin, decreasing to about 2% at 5 GeV. For all 
Pt values above 5 GeV, the correction is < 5%. Systematic uncertainties arising from the 



6 





- 






1 


- 








1 






T 


2 


><10 




sej 


3 












4 












5 


10 










• D0 1994-1996 




b-space (Ladinsky-Yuan) 

CTEQ4Mg,=0.Il GeV' 
gf=0.S8 GeV gj=-l.S GeV' 




— Fixed-order (0(o^)) 

MRSA 



^ 2 







50 



100 150 
Pj<GeV) 



200 



250 



•b-space (Ladinsky-Yuan) 

(x'/d.o.f.=23/26) 

Fixed-Order (0(a;)) 



500 



1 10 10 

P-,(GeV) 

FIG. 2. A comparison of the data to the resummed and fixed-order {0{a1)) calculations. Also 
shown are the fractional differences between the data and the resummed and fixed-order calcula- 
tions. The uncertainties shown include both statistical and systematic uncertainties (other than 
an overall normalization uncertainty due to the luminosity uncertainty). 



choice of ansatz function are evaluated by varying g2 within ±1 standard deviation of the 
best-fit values. Additional uncertainties are evaluated by varying the detector resolutions by 
±1 standard deviation from the nominal values. The effect of these variations is negligible 
relative to the other uncertainties in the measurement. 

Table | shows the values of da{Z e~^e~)/dpT- The uncertainties on the data points 
include statistical and systematic contributions. An additional normalization uncertainty of 
±4.4% arises from the uncertainty on the integrated luminosity [0] that is not included in 
any of the plots nor in the table, but must be taken into account in any fits involving an 
absolute normalization. 

Figure Q shows the final differential cross section, corrected for the D0 detector resolu- 
tions, compared to the fixed-order calculation and the resummation calculation with three 
different parameterizations of the non-perturbative region using published values of the 
non-perturbative parameters. Also shown are the fractional differences of the data from the 
considered resummation predictions. The data are normalized to the measured Z — * e^e~ 
cross section (221 pb |T^) and the predictions are absolutely normalized. We observe the 
best agreement with the Ladinsky-Yuan parameters for the 6-space formalism; however, we 
expect that fits to the data using the Davies- Weber-Stirling (6-space) or Ellis- Veseli (pr- 
space) parameterizations of the non-perturbative functions could describe the data similarly 
well. 

Figure ^ shows the measured differential cross section compared to the fixed-order cal- 
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culation and the resummation calculation using the Ladinsky-Yuan parameterization. We 
observe strong disagreement between the data and the fixed-order prediction in the shape for 
all but the highest values of pt- We attribute this to the divergence of the next-to- leading- 
order calculation at px— 0, and a significant enhancement of the cross section relative to 
the prediction at moderate values of pr- This disagreement confirms the presence of contri- 
butions from soft gluon emission, which are accounted for in the resummation formalisms. 

In summary, we have measured the inclusive difi^crcntial cross section of the Z boson as 
a function of its transverse momentum. With the enhanced precision of this measurement 
over those previous, we can probe non-perturbative, resummation, and fixed-order QCD 
effects. We observe good agreement between the 6-space resummation calculation using the 
published values of the non-perturbative parameters from Ladinsky-Yuan and the measure- 
ment for all values of pr- Using their parameterization for the non-perturbative region, we 
obtain ^2 = 0.59 ± 0.06 GeV^. 

We thank the Fermilab and collaborating institution staffs for contributions to this work 
and acknowledge support from the Department of Energy and National Science Foundation 
(USA), Commissariat a L'Energie Atomique (France), Ministry for Science and Technol- 
ogy and Ministry for Atomic Energy (Russia), CAPES and CNPq (Brazil), Departments 
of Atomic Energy and Science and Education (India), Colciencias (Colombia), CONACyT 
(Mexico), Ministry of Education and KOSEF (Korea), and CONICET and UBACyT (Ar- 
gentina) . 



8 



Pt range 


nominal pT 


number 


aa{Z ^ e^e )/apT 


(GeV) 


value (GeV) 


of events 


(pb/GeV) 


0-1 


0.6 


156 


6.04±0.53 


1-2 


1.5 


424 


16.2±0.96 


2-3 


2.5 


559 


20.4±1.1 


3-4 


3.5 


572 


19.7±1.1 


4-5 


4.5 


501 


16.2±0.92 


5-6 


5.5 


473 


15.0±0.87 


6-7 


6.5 


440 


14.1±0.84 


7-8 


7.5 


346 


ll.lzb0.73 


8-9 


8.5 


312 


10.0±0.69 


9-10 


9.5 


285 


9.29±0.67 


10-12 


11.0 


439 


7.25±0.54 


12-14 


13.0 


326 


5.45±0.44 


14-16 


15.0 


258 


4.45±0.39 


16-18 


17.0 


203 


3.54±0.33 


18-20 


19.0 


181 


3.21±0.31 


20-25 


22.3 


287 


2.06±0.18 




97 3 


1 74 




30-35 


32.3 


124 


0.962±0.11 


35-40 


37.4 


104 


0.840±0.10 


40-50 


44.5 


92 


0.373±0.045 


50-60 


54.5 


61 


0.251±0.036 


60-70 


64.5 


40 


0.163±0.027 


70-85 


76.6 


20 


0.053±0.012 


85-100 


91.7 


13 


0.034±0.009 


100-200 


135 


15 


0.0050±0.0013 


200 300 


228 


2 


oooi+"-"""^ 



TABLE I. Summary of the results of the measurement of the pT distribution of the Z boson. 
The range of pT corresponds to the intervals used for binning the data. The nominal pT corre- 
sponds to the value of pT used to plot the data and was obtained from theory. The quantity 
da{Z — e^e~)/dpT corresponds to the differential cross section in each bin of px for Z — e+e" 
production. The uncertainty on the differential cross section includes both systematic and statis- 
tical uncertatinties, but does not include overall normalization uncertainty due to the luminosity 
of ±4.4%. 
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